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Abstract. We report for the first time a method-independent geometrical expression for the 
angular resolution of an arbitrary network of interferometric gravitational wave (GW) detectors 
when the arrival-time of a GW is unknown. We discuss the implications of our results on how 
to improve angular resolutions of a GW network and on improvements of localization methods. 
An example of an improvement to the null-stream localization method for GWs of unknown 
waveforms is demonstrated. 



1. Introduction 

Several types of astrophysical sources are expected to be detectable both in gravitational waves 
(GWs) and in electromagnetic waves. Coincidence detections of these sources are of significant 
astronomical interest [HE]- A clear understanding of the angular resolution of an array of 
multiple GW detectors is vital to localizations of GW sources and to coincident detections. 

A standard approach to measure how well we can localize a source is to calculate the Fisher 
information matrix where method-independent lower-bounds on statistical errors of estimated 
parameters can be obtained. Numerical results have been calculated by many authors for angular 
resolutions of both the ground-based and the future space GW detector [H HI El [6] . Explicit 
analytical expressions for the angular resolution of a network of GW detectors have been rare 
in the literature. We found two approximate formulae for a 3-detector network which are 
summarized in [7] . One is from private communication of Thorne (cited in Ref. p] ) . The other 
is based on normalized numerical results for a 3-detector network for detections of GWs from 
neutron star-neutron star coalescence using the coherent approach [5j . A general expression for 
an arbitrary network of GW detectors have not been seen. 

Localization of GW sources of unknown waveforms can obtained by the so-called "null- 
stream" method [SI E]- GWs are known to have only two polarizations. The response of an 
interferometric GW detector is a linear combination of the two wave polarizations. Therefore if 
we have data from more than two detectors, we can linearly combine the data to cancel out the 
GW signal. The resulting data streams are called the "null-streams" as they have null-responses 



to GW signals. Localization of a GW source can be achieved by searching for sky directions 
where the constructed null-stream is statistically "null" [U [9] . There are also semi- null streams 
|10j where in linearly combined data, signals are not exactly canceled out but are significantly 
reduced. We propose that localization be further improved by including information from these 
semi-null streams. 

In this report, we summarize results of our on-going research work concerning (1) the angular 
resolution of an arbitrary network of interferometric GW detectors and (2) localization 
methods for GWs of unknown waveforms [12j. An explicit geometrical expression for the angular 
resolution of an arbitrary network of GW detectors is presented for the first time. An improved 
localization method using null-streams combined with semi-null streams is demonstrated and 
compared to that of a straightforward null-stream-only method. 



2. Mathematical Preliminaries 

Suppose we have a network of gravitational-wave detectors, each with spatial size much 
shorter than the GW wavelength, the observed strain of an incoming GW by the individual 
detector / is then a linear combination of the two wave polarizations in the transverse traceless 
(TT) gauge, 

d/(to + r/ + t) = /+/i+(t) + ///ix(t), 0<t<r, (1) 

where to is the arrival time of the wave at the coordinate origin, r/ is the wave travel time 
from the origin to the I-th detector, T is the signal duration, t G [0, T] is the time label of the 
wave. The quantities and are the detector's antenna beam pattern functions [13] for the 
two wave polarizations hx)- They depend on the relative orientation between the detector 
configuration and the frame in which the wave polarizations are defined (which is in turn related 
to the propagation direction n). 

If we assume signal duration to be short enough such that the motion of the detector array 
is unimportant, then in the frequency domain, and in matrix notation, we can write time-delay- 
shifted responses of all detectors as 



d{n) = Ah{n), 

where is the angular frequency. The antenna pattern A is an x 2 constant matrix. 
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and h(r2) is a 2-dimensional vector function. 
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We denote dj as the data from the I-th GW detector and the corresponding noise spectral 
density is 5/, we define a whitened data set of 
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Note that d{0,) corresponds to the whitened data set at each frequency. Correspondingly, we 
denote A as a Nd x 2 response matrix weighted by noise. 
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so that we have d = Ah. For simphcity, we keep the fi-dependence in the notation only when 
it is necessary for clarity. 



3. Geometrical Expression of Angular Resolution 

The angular resolutions are calculated by applying the Fisher information matrix to obtain 
method-independent lower limits on the statistical errors in estimating the direction of a GW 
source. The limits are for unbiased estimators and Gaussian noise (for cautions in using these 
limits, see [14]). The covariance matrix of the "best estimated" angular direction of a GW 
source can be obtained from the corresponding sub-matrix of the inverse of the Fisher matrix 
for all unknown parameters. We show in [11] that in case the initial arriving time to of the 
wave is unknown, the lower bounds of one-sigma error area of angular parameters estimated 
using data from an arbitrary network of GW detectors can be written in a compact geometrical 
form. Here we only summarize the result without showing derivations. We present also only 
cases for short-duration GWs where antenna beam patterns of GW detectors are treated as 
constant. Similar expressions for continuous GWs for ground-based detectors and for the space 
GW detector LISA can be found in ^llj . 

We have defined the error solid angle to be twice the area of the 1-cj error ellipse (measured 
in srad) in angular parameters of 9 (latitude-like) and <f> (longitude-like). 



An = 27r|cos0| \/(A^2^(A02) - (A0A(^)2 , (7) 
we have found that for an arbitrary network of GW detectors, 

'Y.XK,LM ^JK^lm\{ykj X yml) ' np 



where r^j is the displacement vector from detector K to detector J. 

For the worst-case scenario where nothing is known about the initial arrival time to or the 
waveform of a GW, we found that 

1 f+oo ^ ^ 

Akj = - dnn^d*jdkPKj for K^J, (9) 

where matrix P = A(A^A)^"'^A^ (A is defined in Eq. [6]). Note that only J ^ K terms contribute 
in Eq. [51 

For the best-case scenario where the GW waveform is known and the only unknowns are the 
initial wave arrival time to and sky directions, we found 

Akj = for J^K (10) 



where we have defined 



0^2/ °°^n'\dj\'. (11) 
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Note that corresponds to the noise-weighted GW energy fiux coupled to the Jth detector. 



3.1. Implication 

Here we note the clear geometrical meaning of \{tckj x ^ml) ■ n| in Eq. [8l which is twice the 
area of the quadrangle formed by the projections of detectors J, K, L and M onto the plane 
orthogonal to the wave propagation direction. We also note that, in the worst-case scenario 



where the waveform is unknown, the angular resolution is inversely proportional to the weighted 
correlation of responses between detectors (Eq. [9]). In the best-case scenario, it is inversely 
proportional to the fractional GW energy flux coupled to each detector (Eq. [101 Eq. [TTI) . 

Our formula is consistent with the known concept that a larger network is advantageous for 
a better angular resolution. For instance, inclusion of the future Australian AIGO detector can 
improve dramatically the angular resolution of the network [IJ. It further indicates that angular 
resolutions can be improved by optimizing values of A/j. For instance, building more detectors 
of correlated response is advantageous for localizing GWs of unknown waveforms. Similarly, 
a specific localization method can be improved to approach the intrinsic angular resolution 
by selecting data contributing significantly to fractional energy fiux (best-case scenario) or to 
correlations of data between detectors (worst-case scenario). In other words, a localization 
method can be improved by proper treatments of data corresponding to weak responses. 



4. Ranking Network Responses by Singular Value Decomposition Method 

In this and the next section, we demonstrate how one can improve the null-stream localization 
method for GWs of unknown waveforms. We show how to apply the singular value decomposition 
(SVD) method [15] to recombine data from a network of GW detectors to form new data streams 
with characterized sensitivity to GWs from a sky direction. Specifically, we use the SVD to 
construction signal streams, generalized null-streams that have null responses to GWs, and 
semi- null streams that have weak responses to signals. The SVD of A yields (see also ^lOj ) 



A = USVt,S 
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where U and V are unitary matrices of dimensions of N^i x N^i and 2x2 respectively at each 
frequency, i.e., UU"^ = I and VV^ = I, si > S2 > are the so-called singular values. Note they 
are all frequency-dependent. 

We then construct new data streams by inserting this decomposition into equation d = Ah, 
and have 

/ 5i(Wh)i \ 

S2(Wh)2 
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It is evident that the first two components of the new data streams contain signal information, 
the last Nd — 2 terms are null-streams as they have zero-response to signal. In general, null- 
streams can be written as 

/ (Utd)3 \ 
(Utd)4 
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Naturally in case S2 = 0, (UM)2-term should also be included as a null stream. We assume 
for now that S2 ^ 0. The sensitivity level of each signal stream to a GW can be ranked by 
its singular value. Suppose Smax = maxj^QSi(O) {i = 1,2), we define (tentatively) semi-null 



streams as signal streams with corresponding singular values much less than Smax- 



SN(rj) 



if Si <C smax i = 1,2. 



(15) 



Although responses in semi-null streams may not be zero, they can be insignificant compared 
to dominating signal streams. 

5. Improved Localization Strategy Using Semi-Null Streams 

It has been demonstrated [U |9] that null-streams N(r2) can be used to localize a source by 
searching through sky directions for minimum statistic of 



It has also been proposed [TO] that semi-null streams can be included to improve the angular 
resolution. One possible new statistic is 



Instead of searching through sky directions for minimum statistic of Eq. [TO] as discussed in [8l |9] , 
we now search for minimum statistic in Eq. [T7] The tricky part is how to set the threshold at 
which semi- null streams are to be included [12] . 

5.1. Numerical Example 

In this section, results from a Monte-Carlo simulation is presented to illustrate how the null- 
stream localization method can be improved by including the semi-null streams. As a proof- 
of-principle example, we have simply included all semi-null streams in Eq. [T7] that satisfied an 
empirical threshold of Si{Q) / Smax < 0.01 (where Si are singular values defined in Eq. [T2]l . 
Localization is then obtained by searching through sky directions for minimum statistic of 
(1) null-stream-only statistic (Eq. [TO]) and (2) semi-null-stream statistic (Eq. [T7]) respectively. 
Results are then compared. 

We have used simulated signal and noise. For the signal, we used a Sin-Gaussian wave form 
of h-^-{t) = hx{t) = /iQ sin(27r/ot) exp(— i^/r^) with polarization angle chosen arbitrarily to be 
"0 = 0, signal duration T = 7 ms, sampling rate = 16 kHz, central frequency /o = 700 Hz and 
r = 2 ms. The arrival time of the GW wavefront at LIGO Livingston (LI) is chosen arbitrarily 
to be at 0.00 hr, March, 18, 2004. The source direction was chosen to be near that of the 
maximum sensitivity of LI (right ascension RA= 85.1235*^ and declination Dec= 30.56*^) at the 
chosen time. We have chosen an optimal network signal-to-noise ratio of SNR= 20. Location 
information of different GW observatories were obtained from [16] and references therein. For 
the noise, we have adopted the designed noise spectral densities for initial LIGOs (at Livingston, 
LI, and at Hanford, HI) [T7] and for GEO [TO] at 500 Hz tuning. The simulated GW signal 
is then injected into a total of 500 sets of randomly generated Gaussian noise. For each of the 
simulated data of noise plus signal, we use the Nelder-Mead method [19] to search through sky 
directions for minimum statistics of Eq. [TO] and Eq. [T7] respectively. All searches start from 
the source direction to shorten the search time which is adequate for the purpose of proof of 
principle. 

The results are shown in Fig. [T] Source directions obtained using the null-stream-only method 
(Eq. I16|) for different noise realizations are plotted in cross symbols. Filled circles are those from 
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Figure 1. Comparison of localization methods. 



the improved localization method where semi- null streams are included (Eq. [T7|) . The error 
ellipse is for data from the improved method at a 63% confidence level assuming a bi-variate 
normal distribution of angular parameters. The star symbol indicate the average direction in the 
improved method. The actual source direction is indicated with a cross. The triangle symbol 
indicates the average direction from the null-stream-only method. We also plot the time-delay 
lines for Ll-G (gray dotted lines) and Ll-Hl (solid lines) at a 0.2 ms interval. 

It is evident that in this particular example, for a direction where sources are most likely to be 
detected by the LHV-network, inclusion of semi-null streams can improve the source localization 
significantly. The scatter of angular directions obtained from the null-stream-only (gray crosses) 
method is much larger (and therefore worse localization) than that when semi-null streams are 
also included (filled circles). Note that our choice of the semi- null stream is not optimal, further 
improvements are expected when optimal search methods are constructed jl2j . 

6. Conclusion 

We have reported for the first time a compact geometrical expression for an arbitrary network 
of GW detectors when the initial arrival time of a GW is unknown. Our results demonstrate 
the known geometrical elements, as well as the role of energy fluxes and correlation between 
responses of different detectors that determine the intrinsic angular resolution of a GW detector 
array. In the second part of this paper, we show an example where localization of a GW source 
can be improved by including semi-null streams which are linear combination of data that 
have weak response to a GW signal than that of null-stream-only method. We show how the 
Singular- Value-Decomposition method, besides providing a vehicle for generalized optimizations 
of detection methods and for construction of generalized null-streams [lOj, can also be used to 
identify semi-null-streams and help improve the GW source direction determination. 
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